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Abstract 

The injector linacs of the JLC/NLC project include the pre- 
linac, the drive linac, the e~ booster, and the e+ booster. 
The first three will be S-band machines, the last one, an L-band 
machine. We have demonstrated that by using detuning alone 
in the accelerator structure design of these linacs we will have 
acceptable tolerances for emittance growth due to both injection 
jitter and structure misalignments, for both the nominal (2.8 ns) 
and alternate (1.4 ns) bunch spacings. For the L-band struc- 
ture (a structure with 27r/3 phase advance) we take a uniform 
distribution in synchronous dipole mode frequencies, with cen- 
tral frequency / = 2.05 GHz and width A^j = 3%. For the 
S-band case our optimized structure ( a 37r/4 structure) has a 
trapezoidal dipole frequency distribution with / = 3.92 GHz, 
A^j = 5.8%, and tilt parameter a = —.2. The central frequency 
and phase advance were chosen to put bunches early in the train 
on the zero crossing of the wake and, at the same time, keep 
the gradient optimized. We have shown that for random man- 
ufacturing errors with rms 5 /im, (equivalent to 10~^ error in 
synchronous frequency), the injection jitter tolerances are still 
acceptable. We have also shown that the structure alignment 
tolerances are loose, and that the cell-to-cell misalignment tol- 
erance is > 40 /im. Note that in this report we have considered 
only the effects of modes in the first dipole passband. 
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1 Introduction 

A major consideration in the design of the accelerator structures in the 
injector linacs of the JLC/NLC|l][^ is to keep the wakefield effects within 
tolerances for both the nominal (2.8 ns) and the alternate (1.4 ns) bunch 
spacings. One important wakefield effect in the injector linacs is likely to be 
multi-bunch beam break-up (BBU). With this effect a jitter in the injection 
conditions of a bunch train, due to the dipole modes of the accelerator 
structures, is amplified in the linac. By the end of the linac bunches in 
the train are driven to large amplitudes and/or the projected emittance of 
the train becomes large, both effects which can hurt machine performance. 
Another important multi-bunch wakefield effect that needs to be considered 
is static emittance growth caused by structure misalignments. 

To minimize the multi-bunch wakefield effects in the injector linacs we 
need to minimize the sum wake in the accelerator structures. The dipole 
wake amplitude of the structures — and therefore also the sum wake amplitude — 
scales as frequency to the -3 power. Therefore, compared to the main (X- 
band) linac, the injector linac wakes tend to be smaller by a factor 1/64 and 
1/512, respectively, for the S- and L-band linacs. We shall see, however, 
that — in the S-band case — this reduction, by itself, is not sufficient. Two 
ways of reducing the sum wake further are to detune the first pass-band 
dipole modes and to damp them. Detuning can be achieved by gradually 
varying the dimensions of the cells in a structure. Weak damping can be 
achieved by letting the fields couple to manifolds running parallel to the 
structure (as is done in the main JLC/NLC linac stronger damping by, 
for example, introducing lossy material in the cells of the structure. 

In the injector linacs the dipole mode frequencies are much lower than in 
the main linac, and the number of dipole mode oscillations between bunches 
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Ud is much smaller (see Table Therefore, significantly reducing the wake 
envelope at one bunch spacing behind the driving bunch by detuning alone 
becomes more difficult. In addition, for a given Q, the effective damping is 
4 (or 8) times less effective than for X-band. 

Table 1: Scaling of the frequency and the wake amplitude for C, S, and L 
bands as compared to X band. Also given are the number of dipole mode 
oscillations between bunches n^, and the damping needed to reduce the wake 
amplitude by 1/e at the position of the second bunch, Q^, for the nominal 
(2.8 ns) and the alternate (1.4 ns) bunch spacings. 



Scaling 


At = 2.8 ns 


At = 1.4 ns 


Band 


Freq. 


Wake 


rid 


Qd 


nd 


Qd 
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1 
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42.0 


132 


21.0 


66 
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1/2 


1/8 


21.0 


66 


10.5 


33 
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1/4 


1/64 


10.5 


33 


5.3 


16 
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1/8 


1/512 


5.3 


16 


2.6 
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In this report our goal is to design the accelerator structures for the 
injector linacs using simple detuning alone, i.e. including no damping, to 
take care of the long-range wakefields. We focus mostly on the S-band in- 
jector linacs. We begin by discussing analytical approaches to estimating 
the effects of BBU and structure misalignments. We then discuss wakefield 
compensation using detuning. We optimize structure dimensions for struc- 
tures with 27r/3 and 37r/4 per cell phase advance, and show that the latter 
is preferable. And finally we obtain tolerances to wakefield effects for all the 
injector linacs using both analytical formulas and numerical tracking. Note 
that in this report we are only concerned with the effects of modes in the 
first dipole passband, which have kick factors much larger than those in the 
higher passbands. The effects of the higher passband modes, however, will 
need to be addressed in the future. 

2 Emittance Growth 
2.1 Beam Break-up (BBU) 

In the case of single-bunch beam break-up in a linac the amplification of 
injection jitter can be characterized by a strength parameter dependent on 
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the longitudinal position within the bunch. When the strength parameter 
is sufficiently small the growth in amplitude at the end of the linac is given 
by the first power of this parameter!^. For the multi-bunch case we can 
derive an analogous strength parameter, one dependent on bunch number 
m. When this strength parameter is sufficiently small we expect that again 
the growth in amplitude at the end of the linac is given by the first power 
of the parameter. (But even when the strength parameter is not sufficiently 
small it can be a useful parameter for characterizing the strength of BBU.) 
For the multi-bunch case the strength parameter becomes (see Appendix A) 

T^ = ^^^^^^giEj/Eo,0 [m = l,...,M], (1) 

with N the single bunch population, L the machine length, (3q the initial 
value of the beta function averaged over a lattice cell, Eq the initial energy, 
Ef the final energy, and M the number of bunches in a train. The sum wake 
Sm is given by 

m—l 

5^ = Ty[(m - i)At] [m = l,...,M], (2) 

1=1 

with W the transverse wakefield and At the time interval between bunches 
in a train. The wakefield, in turn, is given by a sum over the dipole modes 
in the accelerator structures: 

Wit) = 2kn sin(27r/„t/c) exp(-7r/„t/Q„) , (3) 

n 

with Nm the number of modes, /«, kn, and Qn are, respectively, the fre- 
quency, the kick factor, and the quality factor of the n^^ mode. The function 
g(x) in Eq. |l| is one depending on the energy gradient and focusing proffie 
in the linac. For acceleration assuming the beta function varies as /3 ~ E^, 

5(^,0 = ^(^f^) [iS^E^]. (4) 

If Tm, for all m, is not large, the linear approximation applies, and this 
parameter directly gives the (normalized) growth in amplitude of bunch m. 
If Tm is not large the projected normalized emittance growth of the bunch 
train becomes (assuming, for simplicity, that, in phase space, the beam 
ellipse is initially upright): 

1/2 



6e 



1 + 



CTyO J 



1 [Tm small], (5) 
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with Ho the initial bunch offset, Tj-msO the rms of the strength parameter 
(the square root of the second moment: the average is not subtracted), 
and (7yo the initial beam size. Note that the quantity Sm/M in the multi- 
bunch case takes the place of the bunch wake (the convolution of the wake 
with the bunch distribution) in the single bunch instability problem. As 
jitter tolerance parameter, rt, we can take that ratio y^/cTyo that yields a 
tolerable emittance growth, 5et- 

2.2 Misalignments 

If the structures in the linac are (statically) misaligned with respect to a 
straight line, the beam at the end of the linac will have an increased projected 
emittance. If we have an ensemble of misaligned linacs then, to first order, 
the distribution in emittance growth at the end of these linacs is given by 
an exponential distribution exp[— 5e/(5e)]/(Je), with|^ |^ 

V (<^e) = — ^ ^-h{Ef/Eo,Q , (6) 

with La the structure length, {xa)rms the rms of the structure misalignments, 
Srms is the rms of the sum wake with respect to the average, Na the number 
of structures; the function h is given by (again assuming P ~ E^): 

Eq. 1^ is valid assuming the so-called betratron term in the equation of motion 
is small compared to the misalignment term. 
We can define a misalignment tolerance by 



Xat — \Xa)rms 

with 6et the tolerance in emittance growth. What is the meaning of x^t? 
For an ensemble of machines, each with a different collection of random 
misalignment errors but with the same rms Xat, then the distribution of 
final emittances will be given by the exponential function with expectation 
value 6ef. Note that if we, for example, want to have 95% confidence to 
achieve this emittance growth, we need to align the machine to a tolerance 
level of Xat/V— In .05 « .58xat- 

^This equation is a slightly generalized form of an equation given in Ref. |5|. 
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Besides the tolerance to structure misalignments, we are also interested 
in the tolerance to cell-to-cell misalignments due to fabrication errors. A 
structure is built as a collection of cups, one for each cell, that are brazed 
together, and there will be some error, small compared to the cell dimen- 
sions, in the straightness of each structure. To generate a wake (for a beam 
on-axis) in a structure with cell-to-cell misalignments we use a perturbation 
approach that assumes that, to first order, the mode frequencies remain 
unchanged (from those in the straight structure), and only new kick fac- 
tors are needed Q (The method is described in more detail in Appendix B). 
Note that for particle tracking through structures with internal misalign- 
ments, contributions from both this (orbit independent) wake force and the 
normal (orbit dependent) wake force need to be included. 



Machine properties for the injector linacs used in this report are given 
in Table ^Q. The rf frequencies of all linacs are sub-harmonics of the 
main linac frequency, 11.424 GHz. The prelinac, e"*" drive linac, e~ booster 
linac all operate at S-band (2.856 GHz), and the e"*" booster linac at L-band 
(1.428 GHz). Note that f3yo and (" are only a rough fitting of the real machine 
/3-function to the dependence ^ . In Table ^ beam properties for the 
injector linacs, for the nominal bunch train configuration (95 bunches spaced 
at At = 2.8 ns), are given. For the alternate configuration (190 bunches 
spaced at At = 1.4 ns) is reduced by l/\/2. 

Table 2: Machine properties of the injector linacs. Given are the initial 
energy Eq, the final energy Ef, the length L, the initial average beta function 
in y, and the approximate scaling parameter of /? with energy (/? ~ £"'). 



Name 


Band 


So[GeV] 


Ef[Ge\] 


L[m\ 


/3?;o[m] 


c 


Prelinac 


S 


1.98 


10.0 


558 


8.6 


1/2 


e"*" Drive 


S 


.08 


6.00 


508 


2.4 


1/2 


Booster 


s 


.08 


2.00 


163 


3.4 


1/4 


Booster 


L 


.25 


2.00 


184 


1.5 


1 



3 Wakefield Compensation 

For effective detuning, one generally requires that the wake amplitude drop 
quickly, in the time interval between the first two bunches, and then remain 
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Table 3: Beam properties in the injector linacs under the nominal bunch 
train configuration (95 bunches spaced at Ai = 2.8 ns). Given are the 
bunch population N, the rms bunch length cj^, the initial energy spread 
agQ, and the nominal normalized emittance in y, Cyn- Note that under the 
alternate bunch train configuration (190 bunches spaced at At = 1.4 ns) 
is reduced by 1/ \/2 . 



Name 




az [mm] 


^50[%] 




Prelinac 


1.20 


0.5 


1. 


3 X 10"^ 


Drive 


1.45 


2.5 


1. 


1 X 10~* 


e~ Booster 


1.45 


2.5 


1. 


1 X 10"* 


Booster 


1.60 


9.0 


3.5 


6 X 10"^ 



low until the tail of the bunch train has passed. In the main (X-band) linac 
of the NLC, Gaussian detuning is used to generate a fast Gaussian fall-off in 
the Wakefield; in particular, at the position of the second bunch the wake is 
reduced by roughly 2 orders of magnitude from its initial value. The short 
time behavior of the wake can be analyzed by the so-called "uncoupled" 
model. According to this model (see, for example, Ref. [^) 

Wit)^Yl 2^«n sin(27r/,„t/c) [t small] , (9) 

n 

where Nc is the number of cells in the structure, and fsn and ksn are, re- 
spectively, the frequency and kick factor at the synchronous point, for a 
periodic structure with dimensions of cell n. Therefore, one can predict the 
short time behavior of the wake without solving for the modes of the sys- 
tem. (In the following we will omit the unwieldly subscript s; whether the 
synchronous or mode parameters are meant will be evident from context.) 
For Gaussian detuning the initial fall-off of the wake is given by 

W{t) PS 2fcsin(27r/t) exp (^-2 [irftasf]^^ [t small], (10) 

with k the average kick factor, / the average (first band) synchronous, dipole 
mode frequency, and asf the sigma parameter in the Gaussian distribution. 
Suppose we want a relative amplitude reduction to 0.05 at the position of 
the second bunch. Considering the alternate (1.4 ns) bunch spacing, and 
taking / = 4.012 GHz (S-band), we find that the required (Jsf = 6.5%. 
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To achieve a smooth Gaussian drop-off of the wake requires that we take 
at least ~ 3asf = 20% as the full-width of our frequency distribution, a 
number which is clearly too large. 

If we limit the total frequency spread to an acceptable ^sf = 10% and 
keep the parameter agf fixed, our Gaussian distribution becomes similar to 
a uniform distribution. For the case of a uniform distribution with full width 
Asf the wake becomes 

sin(27r/t) ^^^(^M^/) [(vr/t/g) small]. (11) 

Again considering S-band with the alternate (1.4 ns) bunch spacing, 
and taking / = 4.012 GHz, we obtain an amplitude reduction to 0.37 at 
the position of the second bunch, which is still too large. If we want to 
substantially reduce the wake further we need to shift the average dipole 



frequency /, so that the term sin(27r/t) in Eq. 11 becomes small, and the 



wake at the second bunch is near a zero crossing. That is, 

— Ti 

fAt=— [n an integer] , (12) 

with At the bunch spacing. With n an even integer the bunch train will 
be near the integer resonance, otherwise it will be near the half-integer res- 



onance. With our parameters fAt = 5.62, and Condition [1^ is achieved 
by changing / by —2% (or by a much larger amount in the positive direc- 
tion). However, / = 4.012 GHz is the average dipole mode frequency for the 
somewhat optimized structure, and a change of —2% results in a net loss of 
7% in accelerating gradient, and, presumably, a 7% increase in the required 
lengths of the S-band injector linacs. One final possibility for reducing the 
wake at one bunch spacing is to introduce heavy damping. But for this 
case, just to reduce the wake at one bunch spacing by 1/e, a quality factor 
of 16 would be needed (see Table ||), and such a quality factor is not easy 
to achieve without a significant loss in fundamental mode shunt impedance. 

The Wakefield for a uniform distribution, as given by Eq. not only 
gives the initial drop-off of the wake, but also the longer term behavior. 
(However, here the mode parameters, not the synchronous parameters, are 
needed. Therefore, to see whether such a mode distribution can actually be 
achieved the circuit equations need to be solved.) We see that for a uniform 
distribution the wakefield resurges to a maximum again, at t = Nc/{fAsf). 
Therefore, A^f must be sufficiently small to avoid this resurgence occurring 
before the end of the bunch train; i.e. it must be significantly less than 
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Nc/{M fAt) (which is about 10% in our case). The envelope of Eq. 11 for 
Asf = 8%, / = 4.012 GHz, and = 114 is shown in Fig. |. 

Another possibihty for pushing the resurgence in the wake to larger t is 
to use two structure types, which can effectively double the number of modes 
available for detuning. This idea has been studied; it has been rejected in 
that it requires extremely tight alignment tolerances between pairs of such 
structures. 
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Figure 1: The wake envelope (normalized to a maximum of 1) for a uniform 
frequency distribution. Shown is Eq. |Tl| with all oscillations removed. The 
average frequency / = 4.012 GHz, the bandwidth Asf = 8%, and the num- 
ber of modes Nc = 114. Note that in the injector linacs the bunch train 
extends to ct = 80 m. 



3.1 2n/3 Phase Advance Per Cell 

Except for the region of the initial drop-off, we need to solve for the eigen- 
modes of the system to know the behavior of the wake or the sum wake for 
a detuned structure. To numerically obtain these modes we use a computer 
program that solves the double-band circuit model described in Ref. 1^. We 
consider structures of the disk-loaded type, with constant period and with 
rounded irises of fixed thickness. The iris radii and cavity radii are adjusted 
to give the correct fundamental mode frequency and the desired dipole mode 
spectrum. Therefore, the dimensions of a particular cell m can be specified 
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by one free parameter, which we take to be the synchronous frequency of the 
first dipole mode pass band, fsm (more precisely, the synchronous frequency 
of the periodic structure with cell dimensions of cell m). The computer pro- 
gram generates 2Nc coupled mode frequencies fn and kick factors kn, with 
Nf. the number of cells in a structure. It assumes the modes are trapped at 
the ends of the structure. Only the modes of the first band (approximately 
the first Nc modes) are found accurately by the two-band model. And since, 
in addition, the strengths of the first band modes are much larger than those 
of the second band (in the S-band case the synchronous mode kick factors 
are larger by a factor ~ 35), we will use only the first band modes to obtain 
the wake, and then the sum wake. 

For our S-band structures, we will consider a uniform frequency distri- 
bution, with a central frequency / chosen so that at one bunch spacing, for 
the alternate configuration (At = 1.4 ns), the wake is very close to a zero 
crossing. The strength of interaction with the modes, given by the kick fac- 
tors k, will be stronger near the downstream end of the structure, where the 
iris radii become smaller. To counteract this asymmetry we will allow the 
top of the frequency distribution to be slanted at an angle, and therefore, 
our distribution becomes trapezoidal in shape. We parameterize this slant 
by 

XfifM) + Xfiflo) ' ^ ' 

where A/ is the synchronous frequency distribution, and fio and fhi rep- 
resent, respectively, the lowest and highest frequencies in the distribution. 
Note that —1 < a < 1. With /, a, and the relative width of the distribu- 
tion Asf, we have 3 parameters that we will vary to reduce the wakefield 
effects — specifically by minimizing on the sum wake — for both bunch train 
configurations. 

Each S-band structure operates at a fundamental mode frequency of 
2.856 GHz, and consists of 114 cells with a cell period of 3.5 cm (where the 
phase advance per cell = 27r/3), an iris thickness of 0.584 cm, and cavity 
radius ~ 4.2 cm. The Q of the modes due to wall losses (copper) ~ 14, 500. 
Given our implementation of the SLED-I pulse compression system Q, to 
optimize the rf efficiency the average (synchronous) dipole mode frequency 
needs to be 4.012 GHz. Fig. |^ shows the dispersion curves of the first 
two dipole bands for representative constant impedance, S-band structures, 
with a varying from 1.30 cm to 2.00 cm. The results of a finite element. 
Maxwell Equations solving program, OMEGA2[10|, are given by the plotting 
symbols. The end points of the curves are used to fix the parameters in the 
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circuit program. The two-band circuit results for these constant impedance 
structures are indicated by the curves in the figure. We note good agreement 
in the first band results and not so good agreement in those of the second 
band. For a detuned structure, to obtain the local circuit parameters, we 
interpolate from these representative dispersion curves. 




'O.O 0.1 0.2 0.3 0.4 0.5 



phase advance per cell / 2tt 

Figure 2: The dispersion curves of the first two dipole bands of represen- 
tative constant impedance S-band structures. The phase advance per cell is 
27r/3. Results are given for iris radii of a = 1.30, 1.51, 1.63, 1.80, 1.92, and 
2.00 cm. The plotting symbols give 0MEGA2 results, the curves those of 
the circuit model. The dashed line is the speed of light curve. 

We have 3 parameters to vary in our input (uncoupled) frequency distri- 
bution: the (relative) shift in average frequency from the nominal 4.012 GHz, 
6f, the (relative) width of the distribution A^j, and the flat-top tilt param- 
eter a. Varying these parameters we calculate SrmsO, Srms, and the peak 
value of \S\, for the coupled results, for both bunch train configurations. 
These parameters serve as indicators, respectively, of emittance growth due 
to BBU (injection jitter), emittance growth due to misalignments, and the 
maximum beam excursion due to BBU. From our numerical simulations we 
find that a fairly optimized case consists of 5f = —2.4%, Asf = 7.5%, and 
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a = -0.20. 

In Fig. m we display, for the optimized case, the frequency distribu- 
tion (a), the kick factors (b), and the envelope of the wake (c). The dashed 
curves in (a) and (b) give the uncoupled (input) values. The plotting sym- 
bols in (c) give \ W\ at the bunch positions for the alternate (1.4 ns) bunch 
train configuration. From (b) we note that there are a few modes, trapped 
near the beginning of the structure, which have kick factors significantly 
larger than the rest. This is a consequence of the fact that, for all cells 
of this structure, the dispersion curves are backward waves. From (c) we 
see that, due to these few strong modes, the wake envelope does not nearly 
reach the low, flat bottom that it does for the idealized, uniform frequency 
distribution (see Fig. We note, however, that the short-range drop-off 
is similar to the idealized form (see Fig. |l|), for about 20 m. In addition 
we note that, by setting the second bunch near the zero crossing, many 
following bunches also have wakes with amplitudes significantly below the 
wake envelope. Finally, in Fig. |^ we present the sum wake for both bunch 
train configurations. For this case, for both bunch train configurations, 
SrmsQ = Srms = -02 MV/nC/m^. Note that if we set 5f back to 0, then, for 
the 1.4 ns bunch spacing option, SrmsO becomes a factor of 20 larger. 

3.2 37r/4 Phase Advance Per Cell 

If we would like to regain some of the 7% in accelerating gradient that we 
lost by shifting /, we can move to a structure where the group velocity at the 
synchronous point is less than for the 27r/3 structure (for the same /). One 
solution is to go to a structure with a 37r/4 synchronous point. Note that in 
such a structure the cell length is 3.94 cm and that there are 102 cells per 
structure. Note also that for the same group velocity for the fundamental 
mode a higher phase advance implies larger values of iris radius a, which will 
also improve the short-range wakefield tolerances. The dispersion curves are 
shown in Fig. ^. Note that, in this case, our distribution will have /o < for 
the cell geometries near the beginning of the structure, /o > fn for the cell 
geometries near the end of the structure, while the synchronous phase is near 
pi phase advance. Consequently modes touching either end of the structure 
will only weakly interact with the beam (see, eg Ref. Q), allowing us to have 
a smoother impedance function, and therefore a more uniformly suppressed 
wakefield envelope. This was not the case for 27r/3 structure, where the 
dispersion curves for all cells have a negative slope (between and vr phase 
advance) (see Fig. |2|); it is also not the case for a Sir/ 6 structure, where the 
slopes would all be positive. 
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Figure 3: The optimized 27r/3 structure: 6f = -2.4%, Asf = 7.5%, and 
a = —0.20. Given are the frequency distribution Xf (a), the mode kick 
factors k (b), and the envelope of the wake W (c). The dashes in (a) and 
(b) give the uncoupled results; the plotting symbols in (c) give \W\ at the 
bunch positions for the alternate (1.4 ns) bunch train configuration. Note 
that we display only the modes of the first dipole band, and the wake due 
to these modes. 
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Figure 4: The sum wake for the optimized 2tt/3 structure, for both the 
nominal and alternate bunch train configurations. In both cases SrmsO = 
Srms = -02 MV/nC/m2. 

5.5 I — I — I — I — I — I — \ — \ — \ — I — I — I — I — I — I — I — I — I — I — \ — \ — \ — \ — I — I — 




0.0 0.1 0.2 0.3 0.4 0.5 



phase advance per cell / Stt 



Figure 5: The dispersion curves of the first two dipole bands of represen- 
tative cells in a 37r/4 structure. Results are given for iris radii of a = 1.33, 
1.48, 1.63, 1.80, and 1.98 cm. The plotting symbols give 0MEGA2 results, 
the curves those of the circuit model. The dashed line is the speed of light 
line. 
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Figure 6: For the 37r/4 structure: dependence of 3 sum wake parameters 
on 6f for the nominal (left frames) and alternate (right frames) bunch train 
configurations. The ordinate units are MV/nC/m^. The optimum, 6f = 
—2.3%, Asf = 5.8%, and a = —0.20, is indicated by the plotting symbol. 
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Figure 7: For the 37r/4 structure: dependence of 3 sum wake parameters 
on A^j, for the nominal (left frames) and alternate (right frames) bunch 
train configurations. The ordinate units are MV/nC/m^. The optimum, 
5f = —2.3%, /S.if = 5.8%, and a = —0.20, is indicated by the plotting 
symbol. 
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Figure 8: For the 37r/4 structure: dependence of 3 sum wake parameters 
on Q, for the nominal (left frames) and alternate (right frames) bunch train 
configurations. The ordinate units are MV/nC/m^. The optimum, 6f = 
—2.3%, Asf = 5.8%, and a = —0.20, is indicated by the plotting symbol. 
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Again optimizing on the sum wake, we find that for a fairly optimized 
case 5f = —2.3%, ^sf = 5.8%, and a = —0.20. The change of the indicators 
15*1, Srms and SrmsO as we deviate from this point, for both bunch train 
configurations, is shown in Figs. |6[^ In Fig. ^ we show the / dependence. 
We see that, for both bunch train configurations the results are very sensitive 
to /. In Fig. ^ we give the Asf dependence. We can clearly see the effect of 
the resurgence in the wake when Asf > 7%. And finally, in Fig. |l^ we give 
the a dependence. We note that the tilt in the distribution helps primarily 
in reducing the sensitivity to BBU for the nominal (2.8 ns) bunch train 
configuration. 

In Fig. ^ we display, for the optimized Sir/ 4 case, the frequency distribu- 
tion (a), the kick factors (b), and the envelope of the wake (c). From (b) we 
note that in this case k{f) is a relatively smooth function, as was expected 
from our earlier discussion. From (c) we see that the wake envelope reaches 
a broader, flatter bottom than for the 27r/3 structure, again as we expected. 
Again we note that many of the earlier bunches have wakes with amplitudes 
significantly below the wake envelope. Finally, in Fig. ^ we show the sum 
wake for both bunch train configurations. The rms of these sum wakes are 
much smaller than for the 27r/3 structure: SrmsO = Srms = -004 MV/nC/m^. 

3.3 Frequency Errors 

How sensitive are our results to manufacturing errors? We will begin to 
explore this question by looking at the dependence of the sum wake on 
errors in the synchronous frequencies of the cells of the structure. Note that 
the synchronous frequency of a cell is not equally sensitive to each of the cell 
dimensions. Basically there are 4 dimensions: the iris radius a, the cavity 
radius b, the iris thickness d, and the period length p. The synchronous 
frequency fs is insensitive to d and p, and for the average S-band cell we 
find that Sfg = —.855b and 6fs = —.156a. Or, a —1 micron change in b 
results in Sfs = 2 x 10~^; a —1 micron change in a results in 6fs = 1 x 10~^. 
As for attainable accuracy, let us assume that each synchronous frequency 
can be obtained to a relative accuracy of 10~^, or to about .5 MHz. 

As for systematic frequency errors we note from Fig. ^ that we are es- 
pecially sensitive to changes in average frequency. For example, to dou- 
ble 

SrmsO from its minimum, requires a relative frequency change of only 
4 X 10~^. If each cell frequency has an accuracy of 10~^, and there are 
about 100 cells, the accuracy in the centroid frequency should be ~ 10~^. 
Therefore, the effect of this type of systematic error should be negligible. 

As for random manufacturing errors, let us distinguish two types: "sys- 
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Figure 9: For the optimized 37r/4 structure, 5f = —2.3%, ^Sf = 5.8%, 
and a = —0.20: the frequency distribution Xf (a), the mode kick factors 
k (b), and the envelope of the wake W (c). The dashes in (a) and (b) give 
the uncoupled results; the plotting symbols in (c) give \W\ at the bunch 
positions for the alternate (1.4 ns) bunch train configuration. Note that we 
display only the modes of the first dipole band, and the wake due to these 
modes. 
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Figure 10: The sum wake for the optimized 37r/4 structure, given for 
both the nominal and alternate bunch train configurations. For both cases 
Srmso = Srms = -004 MV/nC/m2. 



tematic random" and "purely random" errors. By "systematic random" we 
mean errors, random in one structure, that are repeated in all structures 
of the prelinac subsystem. "Purely random" errors are, in addition, ran- 
dom from structure to structure. In Fig. 11 we give the resulting SrmsO and 
Srms, for both bunch train configurations, when a random error component 
is added to the (input) synchronous frequencies of the optimal distribution. 
With a frequency spacing of ~ 8 x 10""^, an rms frequency error of 10~^ is a 
relatively small perturbation. We see that for the alternate (1.4 ns) bunch 
spacing the effect of such a perturbation is indeed very small, whereas for 
the nominal (2.8 ns) bunch spacing the effect is large. The reason is that 
with the 1.4 ns bunch spacing the beam sits near a half- integer resonance, 
whereas for the 2.8 ns spacing it sits near the integer resonance. (Resonant 
multi-bunch wakefield effects are discussed in Appendix C.) Note, however, 
that if we consider the case of "purely random" machining errors, with a 
relative accuracy in synchronous frequencies of 10~^, and considering we 
have Nstruc = 140, 127, 41 structures in, respectively, the prelinac, the e"*" 
drive linac, and the e~ booster, then, with a 1/ ^/W^truc reduction in sensi- 
tivity, the appropriate abscissas in the figure become .8, .9, and 1.6 x 10~^. 
At these points, for the 2.8 ns spacing, we see that SrmsO is only a factor 
2±1, 2±1, 3±2 times larger than the zero error result. Finally, as for 
the "systematic random" errors, it is difficult to judge how large they might 
be in the real structure; however, they are likely an order of magnitude less 
than the purely random errors, and should therefore not yield a sum wake 
much larger than that due to the purely random manufacturing errors. 



20 



If we make a weak damping, approximate calculation, by redoing the 
calculation but now with Q = 1000, we find no appreciable effect on the 
resonance behavior for the At = 2.8 ns case with frequency errors. For 
strong damping, taking Q = 100, however, we do find a suppression of the 
resonance effect. 
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Figure 11: The optimal 37r/4 structure: the effect of random frequency 
errors. Shown are the relative (synchronous) frequency error, aj^err vs SrmsO- 
The dashed curves give Srms- Each plotting symbol, with its error bars, 
represents 400 different seeds. 



3.4 e+ Booster 

For the booster (the L-band machine) each accelerator structure consists 
of 72 cells and the synchronous phase advance is taken to be 27r/3. The 
synchronous dipole mode distribution is taken to be uniform, with / = 

2.05 GHz and Asj = 3%. We take Q = 20,000. Note that in this case 
/At = 2.87; for the second bunch to sit on the zero crossing of the wakefield 
would require a shift in frequency of +4.5% or —13%. For the L-band 
structure not much can be gained by changing the frequency spectrum. We 
do gain, however, a factor of 8 reduction in wake in going from S- to L-band, 
which, as we shall see, suffices. 
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4 Tolerances 



For our designed structures we perform the tolerance calculations presented 
in Section 2. To check the analytical estimates, tracking, using the computer 
program LIAR||ll|, was also performed. The numerical simulations were 
simplified in that one macro-particle was used to represent each bunch, and 
the bunch train was taken to be mono-energetic. The analytical BBU results 
and those of LIAR, at the end of the four injector linacs, are compared in 
Table ^. Results are given for bunch spacings of 1.4 and 2.8 ns. Under 
the heading "Analytical" are given the rms of the sum wake Srms > the rms 
of the growth factor Trms, the maximum (within the bunch train) of the 
growth factor T, and the tolerance for 5^t = 10%, i.e. that ratio yo/cTyo 
that results in 10% emittance growth. Under the heading "Numerical" we 
give the LIAR results: the maximum (within the bunch train) growth in 
normalized phase space and the tolerance rt for 10% emittance growth, 
both referenced to the centroid of the first bunch in the train. 

Table 4: Beam break-up calculations for the two bunch spacings. Given 
are the rms of the sum wake SrmsO (in units of MV/nC/m^), the rms and 
the peak of the strength parameter, TrmsO and T, respectively, and the 
analytically obtained tolerance rj for d^t = 10%, i.e. that ratio yo/fyO that 
results in 10% emittance growth. Also shown are LIAR results: the peak 
blow-up in normalized phase space, ^, and the tolerance rt. 



At 


Name 


Analytical 


Numerical 


SrmsO 




T 


n 




n 


2.8ns 


Prelinac 


.004 


.007 


.025 


70 


.020 


85 


e"*" Drive 


.004 


.017 


.066 


25 


.047 


35 


e~ Booster 


.004 


.009 


.036 


50 


.026 


65 


e"*" Booster 


.12 


.119 


.227 


3.8 


.153 


5.5 


1.4ns 


Prelinac 


.004 


.004 


.019 


115 


.015 


140 


e~^ Drive 


.004 


.011 


.049 


45 


.035 


60 


Booster 


.004 


.006 


.027 


80 


.019 


105 


Booster 


.30 


.205 


.379 


2.2 


.257 


3.0 



We note from Table ^ that for the S-band machines, TmsO and are 
both small compared to 1, and that the injection jitter tolerance for 10% 
emittance growth is very large. For the L-band machine, the e"*" booster, the 
tolerances are tighter but still acceptable. We note also that the analytical 
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Table 5: Effect of structure misalignments for the two bunch spacings. Given 
are the rms of the sum wake Srms (in units of MV/nC/m^) and the tolerance 
for structure misalignments for a 10% emittance growth, x^t, both obtained 
analytically (Eqs. IJD and by LIAR. 



At 


Name 


Analytical 


Numerical 


Srms 


Xat[mm] 


Xaj[mm] 


2.8ns 


Prelinac 


.004 


2.9 


3.2 


e"*" Drive 


.004 


100. 


120. 


e~ Booster 


.004 


140. 


170. 


Booster 


.022 


590. 




1.4ns 


Prelinac 


.004 


4.6 


4.8 


Drive 


.004 


150. 


180. 


e~ Booster 


.004 


210. 


260. 


Booster 


.040 


450. 





T agrees well with the numerical ^, as do the two versions of r^. 

In Table ^ we present misalignment results. Given are Srms and the 
tolerance for structure misalignments, Xat, both as given analytically and 
by LIAR. As discussed before, the meaning of Xat is the rms misalignment 
that (for an ensemble of machines) results, on average, in a final emittance 
growth equal to a tolerance 6et, which in this case we set to 10%. From 
Table | we see that the analytical and numerical results agree well, and that 
the misalignment tolerances are all very loose. The tightest tolerance is for 
the prelinac with nominal (2.8 ns) bunch spacing, where the tolerance is still 
an acceptable 3 mm. 

The effect of machining errors will tighten these tolerances for the S- 
band machines with the nominal (2.8 ns) bunch spacing, due to the beam 
being near the integer resonance. If machining adds a purely random error 
component that is equivalent to 10~^ frequency error, we saw earlier that 
(for the 2.8 ns bunch spacing case only) this will tighten the injection jitter 
tolerances by about a factor 2 it 1 for the prelinac and e"*" drive linac, and 
about a factor of 3 it 2 for the e~ booster. But even with this, the toler- 
ances are still very loose. The misalignment tolerances are affected less by 
machining errors. The prelinac tolerance, with 2.8 ns bunch spacing, will 
become ~ 2ibl mm; for 95% confidence in achieving 5^ = 10%, the tolerance 
becomes ~ 1 it .5 mm. 
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Finally, what is the random, cell-to-cell misalignment tolerance? Per- 
forming the perturbation calculation described in Appendix B, and cal- 
culating for 1000 different random structures, we find that Srms = -27 it 
.12 MV/nC/m^ for Ai = 2.8 ns, and Srms = -032 ± .003 MV/nC/m^ for 
At = 1.4 ns. We again see the effect of the integer resonance on the 2.8 ns 
option result. (To verify that this is the case, we performed one run but with 
the bunch spacing changed so that the beam sits near the next half-integer 
resonance (11.5); the result was that Srms dropped by a factor of 6.) For 
the prelinac the cell-to-cell misalignment tolerance becomes 40 /xm for the 
nominal (2.8 ns) bunch configuration and 600 /xm for the alternate (1.4 ns) 
configuration. 

5 Conclusion 

We have demonstrated that by using detuning alone, the four injector linacs 
can be built to sufficiently suppress the multi-bunch wakcficld effects, for 
both the nominal (2.8 ns) and alternate (1.4 ns) bunch spacings. We have 
studied the sensitivity to multi-bunch beam break-up (BBU) and to struc- 
ture misalignments through analytical estimates and numerical tracking, 
and shown that the tolerances to injection jitter, in the former case, and to 
structure misalignments, in the latter case, are not difficult to achieve. We 
have also studied the effect of manufacturing errors on these tolerances, and 
have shown that if the errors are purely random, with an equivalent rms fre- 
quency error of 10^^, then the other tolerances arc still acceptable. Finally, 
we have shown that the cell-to-cell misalignment tolerance is > 40 //m. 

For the L-band machine — the e"*" Booster — we have shown that a uni- 
form detuning of the dipole modes, with central frequency / = 2.05 GHz 
and a total frequency spread A^j = 3%, suffices. For the S-band linacs — 
the Prelinac, the Drive Linac, and the e~ Booster — we have shown that 
the 1.4 ns bunch spacing option forces us to reduce the central frequency 
by 2.3% from the nominal 4.012 GHz. Doing this we lose 7% in effective 
gradient, which, however, can be regained by increasing the phase advance 
per cell from 27r/3 to 37r/4. Our final, optimized distribution is trapezoidal 
in shape with / = 3.920 GHz, A^y^ = 5.8%, and tilt parameter a = —0.2. 

We have demonstrated in this report that the integer resonance, which 
we cannot avoid given the two bunch train alternatives, can make us more 
sensitive to manufacturing errors. Also, we have shown that the analytical, 
single-bunch beam break-up theory, when slightly modified, can be useful in 
predicting the behavior of multi-bunch beam break-up also. Given the rather 
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loose tolerances demonstrated here makes us think that the S-band machines 
can be replaced with C-band ones that still have reasonable tolerances, an 
option which may result in savings in cost, though this needs further study. 
Finally, we should reiterate that in this report we were concerned with the 
effects of the modes in the first dipole passband only. With the wakefield 
of the first band modes greatly suppressed by detuning, the effects of the 
higher bands may no longer be insignificant. This problem will need to be 
addressed in the future. 
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Appendix A: 

Analytical Formula for Weak Multi-Bunch BBU 

In Ref. an analytical formula for single-bunch beam break-up in a smooth 
focusing linac, for the case without energy spread in the beam, is derived, 
the so-called Chao-Richter-Yao (CRY) model for beam break-up. Suppose 
the beam is initially offset from the accelerator axis. The beam break- 
up downstream is characterized by a strength parameter T{t,s), where t 
represents position within the bunch, and s position along the linac. When 
T{t, s) is small compared to 1, the growth in betatron amplitude in the linac 
is proportional to this parameter. When applied to the special case of a 
uniform longitudinal charge distribution, and a linearly growing wakefield, 
the result of the calculation becomes especially simple. In this case the 
growth in orbit amplitude is given as an asymptotic power series in T(i, s), 
and the series can be summed to give a closed form, asymptotic solution for 
single-bunch BBU. The derivation of an analytic formula for multi-bunch 
BBU is almost a trivial modification of the CRY formalism. We will here 
reproduce the important features of the single-bunch derivation of Ref. |^] 
(with slightly modified notation), and then show how it can be modified to 
obtain a result applicable to multi-bunch BBU. Note that we are interested 
in estimating the effect of relatively weak multi-bunch BBU, caused by the 
somewhat complicated wakefields of detuned structures. The more studied 
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multi-bunch BBU problem, i.e. the effect on a bunch train of a single 
strong mode, the so-called "cumulative beam break-up instability" (see, e.g. 
Ref. ||l^), is a somewhat different problem to which our results are not 
meant to apply. 

Let us consider the case of single-bunch beam break-up, where a beam 
is initially offset by distance yo in a linac with acceleration and smooth 
focusing. We assume that there is no energy spread within the beam. The 
equation of motion is 



1 d 
E{s)ds 



E{s) 



dy{t, 



ds 



^'-m = W)' dt'y{t\s)Ht')W{t-t') 



(Al) 

with y[t, s) the bunch offset, a function of position within the bunch t, and 
position along the linac s; with E the beam energy, [l//3(s)] the betatron 
wave number, eNt the total bunch charge, \t{t) the longitudinal charge 
distribution, and W{t) the short-range dipole wakefield. Our convention is 
that negative values of t are toward the front of the bunch. Let us, for the 
moment, limit ourselves to the problem of no acceleration and [5 a constant. 
A. Chao in Ref. Q expands the solution to the equation of motion for this 
problem in a perturbation series 

oo 

y(t,s) = ^y(")(i,5) , (A2) 

n=0 

with the first term given by free betatron oscillation [y^^^ = yocos(s//3)]. He 
then shows that the solution for the higher terms at position s = L, after 
many betatron oscillations, is given by 

.^"^(^,^)-^f^^Vw^/^ , (A3) 

with 



n! V 2E 



Rn{t) = [ dtiX{ti)W{t-ti) r dt2X{t2)W{ti-t2) 
J —oo J —oo 

rtn-1 

dtn X{tn)W{tn-l " tn) , (A4) 



and Ro{z) = 1. An observable y is meant to be the real part of Eq. |A2 . 
The effects of adiabatic acceleration, i.e. sufficiently slow acceleration so 
that the energy doubling distance is large compared to the betatron wave 
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length, and (3 not constant, can be added by simply replacing {(3/E) in 



Eq. A3 by {(5/E), where angle brackets indicate averaging along the linac 
from s = to s = L.^ For example, if the lattice is such that [3 ~ E'^ 
then {13/E) = {po/Eo)g{Ef/Eo,C), where subscripts "0" and "/" signify, 
respectively, initial and final parameters, and 

9i^X) = U^^) [f3-E^]. (A5) 



C V ^ - 1 

Chao then shows that for certain simple combinations of bunch shape 
and wake function shape the integrals in Eq. [A4| can be performed analyt- 
ically, and the result becomes an asymptotic series in powers of a strength 
parameter. For example, for the case of a uniform charge distribution of 
length i (with the front of the bunch at t = 0), and a wake that varies as 
W = W't, the strength parameter is 

2E^e 5(^//^o,C) • (A6) 

If T is small compared to 1, the growth is well approximated by T. If T 
is large, the sum over all terms can be performed to give a closed form, 
asymptotic expression. 

For multi-bunch BBU we are mainly concerned with the interaction of 
the different bunches in the train, and will ignore wakefield forces within 
bunches. The derivation is nearly identical to that for the single-bunch 



BBU. However, in the equation of motion, Eq. Al, the independent variable 
t is no longer a continuous variable, but rather t takes on discrete values 
tm = mAt, where m is a bunch index and At is the bunch spacing. Also, W 
now represents the long-range wakefield. Let us assume that there are M, 
equally populated bunches in a train; i.e. Nt = MN, with the particles 
per bunch. The solution is again expanded in a perturbation series. In the 



solution, Eq. A3, the Rn{t), which are smooth functions of t, are replaced 

by 

^ m— 1 



M 

.7 = 1 



(with 7^° = 1), which is a function of a discrete parameter, the bunch index 
m. Note that TZm'' = Sm/M, with Sm the sum wake. 



Note that the terms Eq. |A3| , related to free betatron oscillation also need 

to be modified in well-known ways to reflect the dependence of (3 on E. ft is the other 
terms, however, which characterize BBU, in which we are interested. 
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Generally the sums in Eq. A? cannot be given in closed form, and there- 
fore a closed, asymptotic expression for multi-bunch BBU cannot be given. 
We can still, however, numerically compute the individual terms equivalent 
to Eq. |A3| for the single bunch case. For example, the first order term in 
amplitude growth is given by 

^n.= '^^^^^9{Ef/E,X) [m = l,...,M]. (A8) 

If this term is small compared to 1 for all m, then BBU is well characterized 
by T. If it is not small, though not extremely large, the next higher terms can 
be computed and their contribution added. For T very large, this approach 
may not be very useful. 

Prom our derivation we see that there is nothing that fundamentally 
distinguishes our BBU solution from a single-bunch BBU solution. If we 
consider again the single-bunch calculation, for the case of a uniform charge 
distribution of length i, we see that we need to perform the integrations 



for Rn in Eq. A4. If we do the integrations numerically, by dividing the 
integrals into discrete steps t„ = (n — l)At and then performing quadrature 
by rectangular rule, we end up with Eq. "Kl with M = i/ ^t. The solution is 



the same as our multi-bunch solution. What distinguishes the multi-bunch 
from the single-bunch problem is that the wakefield for the multi-bunch case 
is not normally monotonic and does not vary smoothly with longitudinal 
position. For such a case it may be more difficult to decide how many terms 
are needed for the sum to converge. 

In Fig. |l2| we give a numerical example: the NLC prelinac with the op- 
timized 37r/4 S-band structure, but with 10~^ systematic frequency errors, 
with the nominal (2.8 ns) bunch spacing (see the main text). The diamonds 
give the first order (a) and the second order (b) perturbation terms. The 
crosses in (a) give the results of a smooth focusing simulation program (tak- 
ing (3 ~ E^^"^), where the free betatron term has been removed. We see that 
the agreement is very good; i.e. the first order term is a good approximation 
to the simulation results. In (b) we note that the next order term is much 
smaller. For this example we find that even if we increase the current by an 
order of magnitude the 1st order term alone remains a good approximation. 
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Figure 12: A numerical example: the NLC prehnac with the optimized 37r/4 
S-band structure, but with 10~^ systematic frequency errors, with the nom- 
inal (2.8 ns) bunch spacing (see the main text). Srms = -005 MV/nC/m^. 
The diamonds give the first order (a) and the second order (b) perturbation 
terms. The crosses in (a) give smooth focusing simulation results with the 
free betatron term removed. 



Appendix B: 

The Wakefield Due to Cell-to-Cell Misalignments 

We assume a structure is composed of many cells that are misaligned trans- 
versely by amounts that are very small compared to the cell dimensions. 
For such a case we assume that the mode frequencies are the same as in the 
ideal structure, and only the mode kick factors are affected. To first order 
we assume that for each mode, the kick factor for the beam on-axis in the 
imperfect structure is the same as for the case with the beam following the 
negative of the misalignment path in the error-free structure. In Fig. ^ we 
sketch a portion of such a misaligned structure (top) and the model used 
for the kick factor calculation (bottom). The sketch is meant to represent a 
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disk-loaded structure that has been built up from a cohection of cups. Note 
that the relative size of the misalignements is exaggerated from what is ex- 
pected, in order to more clearly show the principle. Given this model, the 
method of calculation of the kick factors can be derived using the so-called 
"Condon Method" [jl3| , (see also |l5|). Note that this application to cell- 
to-cell misalignments in an accelerator structure is presented in Ref. [^. 
The results of this perturbation method have been shown to be consistent 
with those using a 3-dimensional scattering matrix analysis]^. We will only 
sketch the derivation below. 
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Figure 13: Sketches of part of a misaligned structure (top) and the model 
used for the kick factor calculation (bottom). Note that the relative size of 
the misalignments here is much exaggerated. 

Consider a closed cavity with perfectly conducting walls. For such a 
cavity the Condon method expands the vector and scalar potentials, in the 
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Coulomb gauge, as a sum over the empty cavity modes. As function of 
position X (x, y, z) and time t the vector potential in the cavity is given as 



A(x,t) = ^(ZA(i)aA(x) , (Bl) 



where 

V2aA + 4aA = , (B2) 

with LOx the frequency of mode A, and a;^ x n = on the metallic surface (n 
is a unit vector normal to the surface). Using the Coulomb gauge implies 
that V • a;^ = 0. The q\ are given by 

(ix + ^lqx = 7^ [ dVj-ax , (B3) 
with the normalization 



eo f 



dVayaix=UxSxy , (B4) 

V 

with j the current density. Note that the integrations are performed over 
the volume of the cavity V. 

The scalar potential is given as 

^y:,t) =J2rx(.t)M^) , (B5) 

X 

where ^ 

V'.^A + %0A = O , (B6) 

with J7a the frequencies associated with and with (^;^ = on the metallic 
surface. The r\ are given by 

^A = ^_^dVp<^A , (B7) 

with p the charge distribution in the cavity. Note that one fundamental 
difference between the behavior of A{x, t) and 'l>(x, t) is that when there are 
no charges in the cavity the vector potential can still oscillate whereas the 
scalar potential must be identically equal to 0. 

Let us consider an ultra-relativistic driving charge Q that passes through 
the cavity parallel to the z axis, and (for simplicity) a test charge following 
at a distance s behind on the same path. Both enter the cavity at position 
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z = and leave at position z = L. The transverse wakefield at the test 
charge is then 

1 ^-^ '■^ 

T,Tr. 2^ 



QLXQ Jq 

QLxq 



dz 



z ~\~ s 

cqx \ — ) y±axz{z) 



z -\- S 

-r\ 1 1 V^<p\{z) 



(B8) 



where the integrals are along the path of the particle trajectory. The param- 
eter xq is a parameter for transverse offset (the transverse wake is usually 
given in units of V/C per longitudinal meter per transverse meter); for a 
cylindrically-symmetric structure it is usually taken to be the offset, from 
the axis, of the driving bunch trajectory. For s > L we can drop the scalar 
potential term (it must be zero when there is no charge in the cavity), and 



the result can be written in the form[15 



W(s) 



E 

A 



2UxujxLxq 



Vx^±Vxe 



>L] 



with 



Vx 



r dzaxz{z)e'^''/'^ 
Jo 



(B9) 



(BIO) 



Note that the arbitrary constants associated with the parameter a;^ in the 
numerator and the denominator of Eq. B9 cancel. Note also that — to the 



same arbitrary constant — | V\ | is the square of the voltage lost by the driving 
particle to mode A and Ux is the energy stored in mode A. 

Consider now the case of a cylindrically-symmetric, multi-cell acceler- 
ating cavity, and let us limit our concern to the effect of the dipole modes 
of such a structure. We will allow the charges to move on an arbitrary, 
zig-zag path in the x — z plane that is close to the axis, and for which the 
slope is everywhere small (so that V_l ~ d/dx). For dipole modes in a 
cylindrically-symmetric, multi-cell accelerator structure, it can shown that 
the synchronous component of axz (the only component that, on average, 
is important) can be written in the form axz = xfx{z) (see e.g. Ref. p6[| ). 



Then Eq. 



19| becomes 



2UxujxLx(] 



i-^xs/c / dz'xiz')fxiz')e'''^^''/^ / dz fx{z)e''^^'/^ 



(Bll) 



>L] . 
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Note that this equation can be written in the form: 

W^{s) = 2A:1 sin ('^ + [s > L] , (B12) 

A ^ 

with k'^ a kind of kick factor and 9x the phase of excitation of mode A. 
Note that in the special case where the particles move parallel to the axis, 
at offset a, k'^ = k\ = c\Vx\'^ /{AUxloxo^L), the normal kick factors for the 
structure, and 9x = 0. For this case it can be shown that Eq. B12| is valid for 



all s > 0[15|. Finally, note that, for the general case, Eq. |B12| can obviously 



not be extrapolated down to s = 0, since it implies that Wx{0) 7^ 0, which we 
believe is nonphysical, implying that a particle can kick itself transversely. 
To obtain the proper equation valid down to s = we would need to include 



the scalar potential term that was dropped in going from Eq. B5 to Eq. B£ . 

Our derivation, presented here, is technically applicable only to struc- 
tures for which all modes are trapped. The modes will be trapped at least 
at the ends of the structure, if the connecting beam tubes have sufficiently 
small radii and the dipole modes do not couple to the fundamental mode 
couplers in the end cells. For detuned structures, like those in the injector 
linacs discussed in this report, most modes are trapped internally within a 
structure, and those that do extend to the ends couple only weakly to the 
beam; for such structures the results here can also be applied, even if the 
conditions on the beam tube radii and the fundamental mode coupler do not 
hold. We believe that even for the damped, detuned structures of the main 
linac of the JLC/NLC, which are similar, though they have manifolds to 
add weak damping to the wakefield, a result very similar to that presented 
here applies. 

To estimate the wakefield associated with very small, random cell-to- 
cell misalignments in accelerator structures we assume that we can use the 
mode eigenfrequencies and eigenvectors of the error-free structure. We ob- 
tain these from the circuit program. Then to find the kick factors we replace 



x{z) in the first integral in Eq. Bll by the zig-zag path representing the neg- 
ative of the cell misalignments, a path we generate using a random number 
generator. The normalization factor xq is set to the rms of the misalign- 
ments. How can we justify using this method for finding the wake at the 
spacing of the bunch positions? For example, for the 37r/4 S-band struc- 
ture, the alternate bunch spacing is only 42 cm whereas the whole structure 



length L = 4.46 m. Therefore, in principle, Eq. Bll is not valid until the 
11**^ bunch spacing. We believe, however, that the scalar potential fields 
will not extend more than one or two cells behind the driving charge (the 
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cell length is 4.375 cm), and therefore this method will be a good approxi- 
mation at all bunch positions behind the driving charge. This belief should 
be tested in the future by repeating the calculation, but now also including 
the contribution from scalar potential terms. 

In Fig. we give a numerical example. Shown, for the optimized 37r/4 
S-band structure (see the main text), are the kick factors and the phases of 
the modes as calculated by the method described in this section. Note that 
Ox is not necessarily small. 
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Figure 14: The kick factors and phases of the modes for a cell-to-cell 
misalignment example. The structure is the optimized 37r/4 S-band struc- 
ture (see the main text). For this example, for the nominal (2.8 ns) bunch 
spacing, Srms = -32 MV/nC/m^. 



Appendix C: 

Resonant Multi-Bunch Wakefield Effects 

It is easy to understand how resonances can arise in a linac with bunch 
trains. Consider the case of the interaction of the beam with one single 
structure mode. The leading bunch enters the structure offset from the axis 
and excites the mode. If the bunch train is sitting on an integer resonance, 
i.e. if fAt = n, with / the mode frequency, At the bunch spacing, and n an 
integer, then when the 2nd bunch arrives it will excite the mode at the same 
phase and also obtain a kick due to the wakefield of the first bunch. The 
mth bunch will also excite the mode in the same phase and obtain (m — 1) 
times the kick from the wakefield that the second bunch experienced (for 



35 



simplicity we assume the mode Q is infinity). On the half- integer resonance, 
i.e. when /At = n + .5, the mth bunch will also receive kicks from the 
Wakefield left by the earlier bunches, but in this case the kicks will alternate 
in direction, and no resonance builds up. For a transverse wakefield effect, 
such as we are interested in, however, this simple description of the resonant 
interaction needs to be modified slightly. For this case the wake varies as 
sin(27r/t), and neither the integer nor the half-integer resonance condition 
will excite any wakefield for the following bunches. In this case resonant 
growth is achieved at a slight deviation from the condition /At = n, as is 
shown below. 

In the following, for simplicity, we will use the "uncoupled" model (which 
is described in Chapter 3 of the main text) to investigate resonant effects in 
the sum wake for a structure with modes with a uniform frequency distribu- 
tion. The point of using the uncoupled model is that it allows us to study 
the effect of an idealized, uniform frequency distribution. As we have seen in 
the main text, an ideal (input) frequency distribution becomes distorted by 
the cell-to-cell coupling of an accelerator structure. As example we will use 
the parameters of a simplified version (all kick factors are equal) of the opti- 
mized 37r/4 S-band structure described in the main text; for bunch structure 
we consider the nominal bunch spacing (At = 2.4 ns). The results for the 
real structure, with coupled modes, will be slightly different yet qualitatively 
the same. Note that we are also aware of a different analysis of resonant 
multi-bunch wakefield effect [ 17 1. 



Consider first the case of a structure with only one dipole mode, with 
frequency /, and a kick factor that we will normalize (for simplicity) to 1/2. 
Suppose there are M bunches in the bunch train. The sum wake at the mth 
bunch is given by 



5«(/At) = ^sin(27r[i-l]/At) 

i=l 

sin (7r[m — 1]/At) sin {nrnfAt) 
~ sin (vr/At) 



(CI) 



As with the nominal (2.8 ns) bunch spacing in the S-band prelinacs, let us, 
for an example, consider M = 95 bunches and the region near the 11th 



harmonic. In Fig. 15 we plot /At vs the sum wake for the Mth (the last) 



bunch, S^^J , near the 11th integer resonance. It can be shown that, if M is 
not small, the largest resonance peaks (the extrema of the curve) are at 

3 

/At -'n'±— [M not small] , (C2) 
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with values ib.72M. Note that at the exact integer and half- integer resonant 
spacings the sum wake is zero. 
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Figure 15: The sum wake at the last bunch in a train vs bunch spacing, 
due to a single mode (Eq. |Cl| ); M = 95 bunches. 

Now let us consider a uniform distribution of mode frequencies. For 
simplicity we will let all the kick factors be equal, and be normalized to 1/2. 
The sum wake, according to the uncoupled model, becomes 



1 



n=l 



/At 1 + 



n 



Nc/2) 



A 



(C3) 



with Nc the number of cells (also the number of modes), / the central fre- 
quency, and Asf the total (relative) width of the frequency distribution. 
As an example, let us consider the optimized Sir /A S-band structure, with 
Nc = 102 and Agf = 5.8%. The sum wake at the last (the Mth) bunch 
position, Sm, is plotted as function of fAt in Fig. Note that the uni- 
form frequency distribution appears to suppress the integer resonance. The 
extrema of the curve (the "horns") that are seen at fAt = 11 it .32 are res- 
onances due to the edges of the frequency distribution, with the condition 
fAt 11/(1 ± '^<5//2)- Note, however, that the sizes of even these spikes 
are small compared to those of the single mode case. 



37 



10.6 10.8 11.0 11.2 11.4 
fAt 

Figure 16: The sum wake at the last bunch in a train vs bunch spacing, due 
to a uniform distribution of mode frequencies (Eq. |C^ ) . The total frequency 
spread Asf = 5.8%, and iV^ = 102. 



Suppose we add frequency errors to our model. We can do this by, in 
each term in the sum of Eq. |C3| , multiplying the frequency by the factor 
(1 + (J/err^'n)) with (5 /err the rms (relative) frequency error and r„ a random 
number with rms 1. Doing this, considering a uniform distribution in fre- 



quency errors with rms 5ferr = 10~ , Fig. 16 becomes Fig. 17. Note that 
this perturbation is small compared to the frequency spacing 5.7 x 10~^, 
so it does not really change the frequency distribution significantly. Never- 
theless, because of resonance-like behavior we can see a large effect on Sm 
throughout the range between the horns of Fig. |l^ (10.68 < fAt < 11.32). 
To model cell-to-cell misalignments, we multiply each term in the sum of 
Eq. by the random factor r„. The results, for a uniform distribution of 
errors with rms 1, are shown in Fig. |l^. Again resonance-like behavior is 
seen throughout the range between the horns of Fig. 1£ . 

We can understand these results in the following manner: Only when 
there are no errors does using a uniform frequency distribution suppress the 
resonance in the region near the integer resonance. But otherwise, using a 
uniform frequency distribution basically only reduces the size of the reso- 
nances, at the expense of extending the range in bunch spacings where they 
can be excited. Instead of being localized in the region near the integer 
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Figure 17: The sum wake at the last bunch in a train vs bunch spacing, 
due to a uniform distribution of frequencies, including frequency errors. The 
total frequency spread A^j = 5.8%, the number of modes N^, = 102, and 
rms relative frequency error is 10~^. 

resonance (/At ~ n), resonance- Hke behavior can now be excited anywhere 
between the limits 

Note that this implies that if A^j > l/(/At), then the resonance-like be- 
havior cannot be avoided no matter what bunch spacing (fractional part) is 
chosen. For example, for the X-band hnac in the NLC, where the total width 
of the dipole frequency distribution (of the dominant first band modes) is 
10%, even for the alternate (1.4 ns) bunch spacing, where the integer part 
of /At is 21, the resonance region cannot be avoided. 
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Figure 18: The sum wake at the last bunch in a train vs bunch spacing, 
due to a uniform distribution of frequencies, including random misahgnment 
errors with rms 1. The total frequency spread A^j = 5.8% and then number 
of modes A^c = 102. 
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